Longitudinal Proximity Effects in Superconducting Transition-Edge Sensors 
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We have found experimentally that the critical current of a square superconducting transition-edge 
sensor (TES) depends exponentially upon the side length L and the square root of the temperature 
T. As a consequence, the effective transition temperature T c of the TES is current- dependent and at 
fixed current scales as 1/L 2 . We also have found that the critical current can show clear Fraunhofer- 
like oscillations in an applied magnetic field, similar to those found in Josephson junctions. The 
observed behavior has a natural theoretical explanation in terms of longitudinal proximity effects if 
the TES is regarded as a weak link between superconducting leads. We have observed the proximity 
effect in these devices over extraordinarily long lengths exceeding 100 pm. 

PACS numbers: 74.25.-q,74.78.Bz,74.25.0p 


A superconductor cooled through its transition tem- 
perature T c while carrying a finite dc bias current 
undergoes an abrupt decrease in electrical resistance 
from its normal-state value Rn to zero. Super- 
conducting transition-edge sensors (TESs) exploit this 
sharp transition; these devices are highly sensitive re- 
sistive thermometers used for precise thermal energy 
measurements. 1 TES microcalorimeters have been devel- 
oped with measured energy resolutions in the X-ray and 
gamma-ray band of A E = 1.8±0.2 eV FWHM at 6 keV, 2 
and A E = 22 eV FWHM at 97 keV, 3 respectively - with 
the latter result at present the largest reported E/AE of 
any non-dispersive photon spectrometer. TESs are suc- 
cessfully used across much of the electromagnetic spec- 
trum, measuring the energy of single-photon absorption 
events from infrared to gamma-ray energies and photon 
fluxes out to the microwave range. 1 Despite these exper- 
imental successes, the dominant physics governing TESs 
biased in the superconducting phase transition remains 
poorly understood. 1 

To achieve high energy resolution it is important to 
control both the TES’s T c and its transition width AT C . 
Because the energy resolution of calorimeters improves 
with decreasing temperature, they are typically designed 
to operate at temperatures around 0.1 K. For a TES, 
this requires a superconductor with T c in that range. 
While there exist a few suitable elemental superconduc- 
tors, the best results have been achieved using proximity- 
coupled, superconductor/normal- metal (S/N) bilayers 2 ’ 3 , 
for which T c is tuned by selection of the thicknesses of 
the S and N layers. 4 

There have been a variety of models 4-8 used to explain 
the noise, T c , and AT C in TES bilayers, all assuming 
spatially uniform devices. Though some have been shown 
to be consistent with certain aspects of particular devices, 
they do not explain measured T c and AT C in S/N bilayer 
TESs generally. 

In this paper we emphasize the importance of a phe- 
nomenon that so far has been neglected in previous theo- 
retical studies of TESs: the longitudinal proximity effect. 


Since the square bilayers at the heart of the TES are con- 
nected at opposite ends to superconducting leads with 
transition temperatures well above the intrinsic transi- 
tion temperature of the bilayers, superconductivity is in- 
duced longitudinally into the bilayers via the proximity 
effect. As we shall explain later, many of the basic prop- 
erties of our TES structures are well described by regard- 
ing them as SS'S or SN'S weak links. 9-11 

In this paper we report the properties of TESs based 
on square ( L x L) electron-beam-deposited Mo/Au bilay- 
ers consisting of 55 nm Mo layers (T c ~ 0.9 K) to which 
210 nm of Au is added. The square side lengths L range 
from 8 /im to 290 /zm, and the normal-state resistance 
per square is Rn = 17. 2 ±0.5 mfl. The bilayers are con- 
nected at opposite ends to Mo/Nb leads having measured 
superconducting transition temperatures of 3.5 and 7.1 
K. 12 Further details on the device fabrication process can 
be found in Ref. 13. 

Our measurements are made in an adiabatic demagne- 
tization refrigerator (ADR) with mu-metal and Nb enclo- 
sures providing magnetic shielding for the TES devices 
and SQUID electronics. The magnetic field normal to 
the TES device plane is controlled by a superconducting 
coil with the field value determined from the coil geome- 
try and current. Measurements of the TES resistance R 
are made by applying a sinusoidal current of frequency 
5-10 Hz and amplitude Rias ~ 50-250 nA, with zero dc 
component, to the TES in parallel with a 0.2 niQ shunt 
resistor ( R s h ). The time-dependent TES current is mea- 
sured with a SQUID feedback circuit with input coil in 
series with TES. When Ri as is less than the TES crit- 
ical current I c , R is zero, and all the ac current flows 
through the TES. However, when R ias > I c and R > 0 
during part of the ac cycle, the TES current becomes 
non-sinusoidal, and its maximum value I becomes less 
than Ibias- The TES resistance R at the TES current I 
is then determined from R = R s h(Rias ~ I)/I. 

The critical current I c is measured, with the ADR held 
at constant temperature, by ramping the dc bias current 
from zero and defining I c as the TES current at the first 
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measured finite resistance ( R ~ 10 /ifi) across the TES. 
Record averaging is used at higher temperatures where 
I c becomes small. 

The solid curves in Fig. 1(a) and (b) show measure- 
ments of the critical current I c over seven decades vs 
temperature T. Note that although we find the in- 
trinsic transition temperature of the Mo/Au bilayer is 
T cw = 170.9±0.1 inK, at very low currents a zero- 
resistance state is measured up to much higher temper- 
atures as the TES size is reduced, three times T cw for L 
= 8 ^m. On the other hand, for the larger TES sizes 
(L=130 and 290 ^m) the critical current / C (T) decreases 
rapidly with T near T cw . The observed I c behavior as 
functions of both T and the length L provides strong ev- 
idence that our TESs behave as weak-link devices. The 
dotted curves in (a) and (b) show calculated values of 
I c using the Ginzburg-Landau theory described below. 
In addition, at appropriately chosen temperatures, the 
critical currents of these devices exhibit Fraunhofer-like 
oscillations as a function of an applied magnetic field, be- 
havior characteristic of Josephson weak links. 10,14-16 See 
Fig. 1(c) for an example. 

Because I c depends upon T and A, the effective tran- 
sition temperature T c of the TES (the temperature at 
which an electrical resistance first appears, i.e. , R ~ 
10 /ifi) is both current-dependent and length-dependent. 
Figure 2(a) exhibits these effects. The points labeled 
T C (J, L) are the effective transition temperatures at five 
different current levels (10 nA to 100 fiA) for the data in 
Fig. 1(a) and (b), showing that T c — T cw for each current 
level scales approximately as 1 / L 2 (solid curve fits) for L 
ranging from 8 to 290 /im. For each L, T c — T cw depends 
upon the current. 

Also shown in Fig. 2(a) are temperatures Tr- q.ir n 
and Tr=o.5R n f° r which the resistances are R = 0.1 Rn 
and 0.5 Rn, respectively, from which we define A Tr = 
Tq. 5 r n — To.iiijy We also define transition widths from 
the I c measurements AT C i = T c ioo nA — T c io n A and 
A T c2 = T cl/M - T cW nA ■ In Fig. 2(b) we show that 
these three measures of the transition width all vary ap- 
proximately as 1/A 2 , shown by the dotted line. It also 
follows that T c — T cw scales linearly with the transition 
width. 

Weak links in various SS'S or SN'S weak-link struc- 
tures have been studied experimentally and theoretically 
by numerous authors. However, here we use a simple ver- 
sion of Ginzburg-Landau (GL) theory 17,18 to explain the 
results shown in Figs. 1 and 2. This theory describes the 
complex order parameter ip(r ), whose absolute square 
\ip(r)\ 2 is the superfluid density in the weak link. We 
employ the substitution ip = ip r fe lJ , where ip r is the 
magnitude of the order parameter at the reference points 
x = ±L/2 adjacent to the leads, / = \ip\/ip r is the nor- 
malized order parameter, and 7 is the phase. At the 
reference points, the local value of ip 2 is inversely pro- 
portional to the square of the local penetration depth X r 
via ip 2 = m/4/j,oe 2 A 2 , and a characteristic reference cur- 
rent density j r can be defined via j r = </o/27TA t oA 8 , where 
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FIG. 1: (Color online) (a) and inset (b) Measured (solid lines 
and markers) and theoretical (dotted curves) critical current 
I c versus temperature T for square TESs with side lengths 
L ranging from 8 to 290 fi m. The bold continuous segments 
at the lowest currents are obtained by record averaging. The 
intrinsic transition temperature of the Mo/Au bilayer weak 
links is T cw = 170.9±0.1 mK (thin vertical lines). For T 
somewhat larger than T cw , I c decays approximately exponen- 
tially with the square root of T—T ow . I c also depends strongly 
upon L, which is particularly noticeable for the smaller de- 
vices. T and L values of the constant current contours (hori- 
zontal dashed lines) are plotted in Fig. 2(a). Inset (c) shows 7 C 
vs applied field for the L = 29 /im device showing Fraunhofer- 
like oscillations, similar to those seen in Josephson junctions, 
providing further evidence that the TES exhibits weak-link 
behavior. 


(po = h/ 2 e. 

Near the center of the weak link, where / becomes 
very small, the local penetration depth A = A r /f be- 
comes very large. Moreover, in a thin film of thick- 
ness d < A, magnetic fields and currents spread out over 
the two-dimensional screening length (or Pearl length 19 ) 
A = 2A 2 /d = 2\ 2 /df 2 . For all of our samples there is a 
range of temperatures T far enough above T cw that at 
the center of the weak link A>L and the current den- 
sity is j = xj x = xl/Wd , where / is the TES current. 
The first and second GL equations given in Refs. 17 and 
18 then depend only upon x and can be written as 

_ r+ <l_L / + f! / » + ? =0 a, 

Sio /v r /v r J 

and 

j = jx/jr = ~ A r / 2 (y' + 2ttA x / ( po), (2) 
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FIG. 2: (a) Measurements of the effective transition temper- 
ature T c at different currents and lengths. Markers T C (I,L ) 
give the effective T c from constant current contours of the 
I C (T,L) data in Fig. 1, with solid curves being 1 / L 2 fits for 
each current level. Markers T(R) give temperatures where 
R = 0.1 Rn and 0.5 Rn- (b) Three different measures of the 
transition width defined as differences between pairs of cor- 
responding points shown in (a), as labeled A Tr, Alii, and 
AT C 2 , showing 1 /L 2 scaling (dotted line). 


where t = T/T cw is the reduced temperature, £(T) = 
Wl t - II 1 / 2 is the temperature-dependent coherence 
length, k is the dimensionless Ginzburg-Landau param- 
eter, A x is the vector potential, and the primes denote 
derivatives with respect to x. 

In this paper we are concerned chiefly with weak-link 
behavior for which /( x) is an even function of x and has 
a minimum in the middle of the weak link, /( 0) = /o, 
where /'( 0) = 0. We can obtain an equation that deter- 
mines how /o depends upon L , t, and j by multiplying 
Eq. (1) by /', integrating the result, multiplying by / 2 , 
and taking the square root, which yields the following 
equation, valid for 0 < x < L/ 2, 

df 2 /dx = P(f 2 ), (3) 


where 


r/G) d p 

JfS 


( 6 ) 


The gauge-invariant phase difference across the weak link 
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The integrals in Eqs. (5), (6), and (7) can be evaluated 
numerically as in Ref. 20 or in terms of elliptic integrals 
as in Refs. 21 and 22. For given values of A r , n, £ w , t, 
and L, the solutions of Eq. (5) reveal that j is a single- 
valued function of /o, starting with the value j = 0 at 
/o = 0, initially increasing linearly with fo, rising to a 
maximum value defined as j c , then returning to zero at 
a larger value of fo- 

When t > 1, the above equations reveal that j((f>) is a 
single-valued function of <j> and has a functional depen- 
dence close to j = j c sin cj>, similar to that of a Josephson 
junction. For 0 < <f> < n, the reduced order parameter 
f(j, x) at x = 0 has its maximum value /oo = /( 0,0) 
when <j) = 0, its minimum value 0 when </> = n, and 
a value between these two limits at the critical current 
when j = j c and <f> ft 7 t/2. 

When T > T cw and L > £(T) = £ w /pt — 1, / « 1 
for a large fraction of the length L 1 and one may omit 
the term proportional to / 3 on the right-hand side of Eq. 
(1) to obtain the reduced order parameter f(j, x). In the 
absence of a current, f(0,x) = f r cosh(x/£)/ cosh(L/2£) 
near the center of the weak link, f 00 = /( 0,0) = 
f r / cosh(L/2£) « 2 f r e~ L l 2( ‘ at the center, and the gauge- 
invariant phase difference across the weak link is 4> = 0. 
The parameter f r , which is of the order of unity, would 
be equal to unity if the linearized GL equation were valid 
over the entire length L of the weak link; the suppression 
of f r below unity occurs because the exact solution for 
/( 0, a;) near x ft ±L/2 is strongly influenced by the term 
(K/A r ) 2 / 3 on the right-hand side of Eq. (1). 

For nonzero current, /o = f(j, 0), the reduced order 
parameter at the center of the weak link, is suppressed 
below /oo, and the gauge-invariant phase difference <fi 
across the weak link obeys sin(</>/2) = fo/foo- The re- 
duced current is given by 


3 = 


•V j -2 • X 4/ 2 A r ( fo\ [[ ( fo \ 2 -L/£ ro\ 

-, mSin4 ,K— «, , 8 ) 


such that the reduced critical current is given for T > T cw 
and any L by the approximation 


(4 ) 

such that fo and f(x) can be obtained from the integrals 


f 1 dp 

Jfs np) 



( 5 ) 


Jc = jc/ Jr = (A r /2£)/oo ft ( 2p r \ r p)e~ L * (9) 

at the maximum of /, where fo = / 00 /V 2 and (j) = n/2. 
From Eq. (9), we may obtain the critical current as 
I c = j c Ld = jrj c Ld. Inferring T cw = 170.9 rnK from 
the experimental data in Fig. 1(b) for L = 290 pin and 
assuming n = \ r /£ w , we obtained / w and A r by fitting 
the experimental I c data for L = 8 pm at 250 rnK and 
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375 mK. The dotted curves in Fig. 1(a) and (b) show I c 
calculated using = 738 nm and X r = 79 nm in Eq. (9) 
and Eq. (5), from which /oo was obtained. 23 

Under conditions for which Eq. (9) is valid, if we define 
the effective transition temperature T c (j) as the temper- 
ature at which the first voltage appears along the length 
of the TES when it carries a reduced current density 
j, we can_ determine T c (j) or t c (j) = T c (j)/T cw by set- 
ting j = j c in Eq. (9) and solving for t c {j), noting that 
£ = £w/Vt — 1- The result is 

(T c — T cw )/T cw = (^/L 2 )ln 2 (2/ 2 A rV / U^l/j^). (10) 

Since the dependence upon t c on the right-hand side is 
very weak, because it appears within the argument of the 
logarithm, Eq. (10) predicts that the current-dependent 
transition temperature of the TES should scale very 
nearly as T c — T cw oc 1/T 2 and that T c should increase 
as the square of the logarithm of the inverse TES cur- 
rent. Similar reasoning leads to the conclusion that both 
AT) and AT 2 scale as 1/T 2 . Scaling of A Tr can be un- 
derstood using a simple model of the resistive transition 
based on the assumption that R = {2xj / L)Rn, where, 
for a given reduced current density j, Xj is the solution 
of j = (A r /20/(0,x,) 2 . 

We conclude that TESs behave as weak links. This 
conclusion is based on our experimental findings that (a) 
the critical current at the first onset of a voltage along the 
length depends exponentially upon the length L and the 


square root of the temperature T, (b) both the current- 
dependent effective transition temperature T c and the 
transition width scale as 1/T 2 , and (c) the TESs show 
clear Fraunhofer oscillations as a function of applied mag- 
netic field, characteristic of Josephson weak links. It fol- 
lows that the strength of superconducting order is not 
uniform over the TES. Our findings have implications on 
TES magnetic field sensitivity, which impacts required 
limits on ambient magnetic field magnitude and fluctu- 
ations in TES applications. Proposed uses of the lon- 
gitudinal proximity effect for TES applications include 
(1) tuning the effective T c of TES arrays by changing 
L in mask design, which could compensate for bilayer 
T cw variability 1 and increase yield, and (2) making small 
TESs consisting of superconducting leads separated by 
normal metal, such as Au with T r „, = 0, avoiding the use 
of S/N bilayers. 
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